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Abstract 

Renormalization of composite fields is employed to suppress the statistical 
noise in lattice gauge calculations. We propose a new action which differs 
from the standard Wilson action by "irrelevant" operators, but suppresses 
the fluctuations of the plaquette. We numerically study the Creutz ratios 
and find a scaling window. The SU(2) mass gap is estimated. We prove that 
the contributions of the "irrelevant" operators to the screening mass decrease 
towards the continuum limit. The results obtained from the action with noise 
suppression are compared with those of the standard Wilson action. 
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1 Introduction 



The present lattice calculations provide the only rigorous approach to low energy 
Yang-Mills theories. After the scaling window had been discovered by Creutz in 
his pioneering work []T|, lattice simulations provided first informations on the ratio 
of the low lying glue-ball masses and the string tension for the SU(2) [0 and the 
SU(3) I, |] gauge group. 

Unfortunately, the limited capacity of the computers put severe constraints on the 
accuracy of "lattice measurements". Firstly, the finite number of lattice links cor- 
respond to a finite physical volume. Nowadays, a physical volume of (1.6 fm)"^ is 
available for reasonable values of the lattice spacing (see e.g. @]). Secondly, the 
finite number of independent configurations which are employed to calculate the 
expectation value of the desired operator implies that the "lattice measurement" is 
contaminated with statistical noise. 

Weisz and Symanzik have shown that the situation corresponding to the finite size 
problem can be significantly improved by using an improved action In the 
numerical simulation, an effective action is used which already contains corrections 
from perturbative radiation. In recent years, much work has been devoted to the 
development of such improved lattice actions, which are often referred to as "perfect" 
lattice actions 0. 

In this paper, we will focus on the noise problem. In order to outline the conceptual 
nature of the noise problem, we briefly review the arguments presented in . Glue- 
ball (screening) masses rrig are extracted from correlation functions, i.e. 

C{t) := mt) $(0)) ^ const, e"'^''* , (1) 

where the brackets indicate an average over independent lattice configurations. The 
statistical error of the desired quantity is measured by the standard deviation 

(<l'(t)$(0)<l'(t)$(0)) - C\t) ^ ($2(0)) . (2) 

From perturbation theory, one knows that composite operators acquire new di- 
vergences implying that the statistical error of the correlation function C{t), i.e. 

($^(0)), diverges in the continuum limit a — 0. The disastrous and fundamental 
problem therefore is that the signal-to-noise-ratio vanishes in the continuum limit. 
In the recent past, two concepts have been established to be important in order 
to improve the signal-to-noise-ratio. Firstly, the choice of a non-local operator 
$(x) in (|l|) (which nevertheless carries the quantum numbers of the state under 
investigation) might result in a composite operator ^"^{x) which is free of ultra- 
violet divergences [smearing [|]). Secondly, informations from the links of a former 
update step is used to enhance the signal {fuzzing [Q). Hybrid algorithms, which 
combine "smearing" and "fuzzing", as well as an estimate of their impact on the 
signal-to-noise-ratio can be found in 0. 
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In this paper, we prose a new method to improve the signal-to-noise-ratio. The 
method is inspired from the renormahzation procedure for composite operators in 
continuum quantum field theory. The basic idea is to add to the Wilson action an 
additional term which vanishes faster than the Wilson action in the continuum limit, 
i.e. we add an "irrelevant operator", but which suppresses the statistical noise. We 
study the efficiency of our method by calculating the SU(2) mass gap employing the 
"old" idea of plaquette-plaquette correlations. 

The paper is organized as follows: in the next section, we briefly review the renormal- 
ization of composite operators in continuum quantum field theory. We then discuss 
the modifications of the Wilson action by "irrelevant" composite fields which yield 
the suppression of the noise. In the third section, we present our numerical results. 
Conclusions are left to the final section. 



2 Renormalization of composite operators 
2.1 In continuum quantum field theory 

For illustration purposes, we here consider the continuum quantum field theory of 
a field (f){x) which is described by the Euclidean partition function 

Z[r]]{g) = J Vct>{x) exp |-5[</)]((?) + j (fxr^{x)<P{x)^ , (3) 

where a regularization is understood in order to make @ well defined. For simplicity, 
we assume that the Euclidean action S contains only one parameter g (e.g. the 
coupling strength). The external source rj linearly couples to the field 0(x) implying 
that functional derivatives of Z[ri\{g) with respect to r] yield connected Green's 
functions 

(T0(xi)...0(x^)) . (4) 

These Green's functions are generically divergent in four space-time dimensions, if 
the regulator is removed. Renormalized Green's functions are obtained from the 
generating functional 

ZR[r]R]{gR) := Z[Z^r]R\{ZggR) (5) 

by performing the functional derivative with respect to 1]^. The Zs are the so-called 
renormalization constants, and tir is the renormalized source which accounts for 
field renormalization, and g^ is the renormalized parameter. In order to guarantee 
that (j^) yields finite Green's functions, we have tacitly assumed that the field theory 
(!]) is multiplicative renormalizable |]10|, i-e. all divergences can be absorbed in the 



renormalization constants Z. 

The crucial observation is that, if we allow for composite field insertions, i.e. if we are 
interested in the limit Xi X2 in the renormalized Green's function, new divergences 
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arise in Zji[ri^]{gji). In order to renormalize these insertions, we generalize 



ZR[riR]{gR) ZR[r]n,jR]{gR) = (6) 
V(f){x) exp <^ -S[(t)]{Zrg) + / d'^x Z^r]R{x)(l){x) + ZjjR{x)(f)^{x) 



The additional divergences due to the composite field can be absorbed in the 

renormalization constant Zj. 

The dependence of the renormalization constants on the regulator is of course not 
known a priori. Perturbation theory usually provides a systematic way to extract 
this dependence [0. In the context of numerical lattice gauge calculations, the 
question arises, how one should choose the bare source j{x) as function of the lattice 
spacing in order to renormalize the composite field insertions and therefore to reduce 
the statistical noise in the continuum limit. We will answer this question in the next 
section. 



2.2 In lattice gauge calculations 

The partition function of SU(2) lattice Yang-Mills theory is defined as a functional 
integral over the link variables U^{x), i.e. 

Z = JVU, exp{-S} , (7) 

where the standard Wilson action is given by 

Sw = E /5 [1 - PA^)] ■ (8) 

P^j^y is the plaquette, which is built from four link variables, i.e. 

P^u{x) = itr{f/^(x)f/,(x + /i)f/t(x + z/)f/t(x)} . (9) 

The functional integral (0) is defined on a lattice with lattice spacing a, which serves 
as the ultraviolet regulator. In the continuum limit (a — >^ 0), the plaquette is 

PA^) = 1 - j^;.^;. + o{a') , (10) 

where F^^ is the usual field strength tensor. In this paper, we will confine ourselves 
to the plaquette-plaquette correlation function (P^,^(x)Pa/3(0)) in order to extract 
the mass gap of the SU(2) lattice theory. From (|^), it is clear that this correlation 
function is plagued by a statistical noise which diverges in the continuum limit. 
From the discussions in the last section, it is now evident that one must add a term 
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J2{x}iiu j{^)P^ui^) "with a suitable choice of j{x) to the action (|]) in order to avoid 
this divergence. We here propose to perform the numerical simulation using the 
action 

s = j: p[i- p,.(x)] + y: j {p,uix)-Af , (11) 

{x}fj,u {x}fj,u 

where j and A are constants. In fact, we will choose A to be the average value of 
the plaquette Pf^^, which will be a function of /? and j. 

Let us study the naive continuum limit of the action 5" (pA]). Using (|l3), a direct 
calculation yields (up to a constant) 

S = [P- 2j (1 - A)] + 0{a') . (12) 

This implies that the action S cannot be distinguished in the naive continuum limit 
from Wilson's action (H) with an effective parameter /3ejj = (3 — 2j{l — A). In 
the quantum continuum limit (/? oo), the average plaquette and effective inverse 
temperature Peff are approximately given by |l[| 

^ = A/z^f-l- (13) 

If j increases less than linearly with (3, the results of the quantum theory using 5* 
should agree with those which are obtained by employing the Wilson action. 
On the other hand, the term in ([Tl|) proportional to j further constrains the plaquette 
to its average value A and therefore suppresses statistical fluctuations around the 
average value of the plaquette. The key point is that the action 5* (0) differs 
from the Wilson action by "irrelevant" operators which are chosen to suppress the 
statistical noise of the plaquette. The prize one has to pay is that the average 
plaquette value must already be known at the beginning of the numerical simulation. 



3 Numerical results 
3.1 Creutz ratios 

We perform our numerical simulations of the quantum theory employing the action 
(PH) on a 10"^ lattice using the heat bath algorithm by Creutz [Q. Our purpose is to 
demonstrate the mechanism of noise suppression proposed in the previous sections, 
rather than to provide new precision measurements. In the latter case, one should 
resort to "improved" actions as well as a larger number of lattice points. 
The first task is to calculate the average plaquette A self-consistently for given values 
for P and j. We apply the following procedure (we leave it to the reader to develop 
his own method): before the lattice has reached its thermo-dynamical equilibrium. 
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Figure 1: The Creutz ratios employing the standard Wilson action (left) and using 
the action S (|l^) with j = 0.5 (right). The dashed line in the right hand picture 
indicates the perturbative scaling behavior of the case j = 0. 



we use the lattice average from the previous heat bath step for A. When equilibrium 
is reached, we calculate the average plaquette taking into account all heat bath steps. 
In a particular step, we assign the actual value of this average plaquette to A. We 
then perform a large number of heat bath steps to obtain an accurate value of A, 
which subsequently enters the numerical calculations of correlation functions, where 
a smaller number of heat bath steps is sufficient. 

The crucial question is whether the scaling limit is reached for finite values of j . In 
order to answer this question, we calculate the Creutz ratios as a function of f3. 
The left picture of figure |l| shows the case j = 0. These are the Creutz ratios which 
one obtains using the standard Wilson action. These results are compared with those 
for j = 0.5 (in the right picture of figure 0). The lines indicate the scaling behavior 
which has been calculated with the help of the perturbative renormalization group (3- 
function. The crucial observation is that the model with non-vanishing source j also 
approaches the scaling behavior (shown by the lines in figure ^ which is predicted 
by the perturbative renormalization group. Also note that the perturbative scaling 
already sets in at smaller values of j3 compared with those of the case j = 0. This 
is precisely what one expects from (|T^) . 

3.2 Noise suppression 

Let us study the efficiency with which the action in ([TT|) suppresses the statistical 
error of the average value of the plaquette. For this purpose we consider the proba- 
bility distribution of finding a particular value of the plaquette P^i, in the interval 

- 5 X 10"^ ^ + 5 X 10^3] . (14) 
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Figure 2: The distribution of the plaquette around its average value A for several 
values of the noise suppression factor j. 



From the numerical point of view, we proceed as follows: we divide the above 
interval in bins of width 10""^, and calculate the lattice average of the plaquette in 
a particular heat bath step. We then count the number of average values which 
correspond to a certain bin. We evaluated 1140 heat bath steps. The numerical 
result is shown in figure ^ for /3 = 2.1 and j = 0.5. One clearly observes that the 
data points are strongly grouped around the corresponding average value A for large 
values of the noise suppression factor j. In addition, one observes a sharp peak at 
Pfiu = A. Whether this peak is an artifact due to corrections to the action of order 
or whether the peak is necessary to reproduce the standard Yang-Mills action 
(0) in the scaling limit /5 — cxd, is not clear to us. 

3.3 The SU(2) mass gap 

In this subsection, we will numerically estimate the SU(2) mass gap from the 
plaquette-plaquette correlation function in order to demonstrate how the action 
(PH) works in practice. The purpose of this subsection is two-fold: Firstly, we want 
to show that the mass gap obtained here is in agreement with the high precision 
measurements which employ high statistics and an improvement of the signal- 
to-noise-ratio using the fuzzing technique. Secondly, we will compare the results at 
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several values of (3 in order to estimate the contribution of the "irrelevant" operators 
off the continuum limit. We are aware of the problem that the overlap of the pla- 



quette with glue-ball wave-function is small [|1^]. For high precision measurements, 
one should therefore employ non-local operators. Performing the noise suppression 
for the case of these operators, however, might be numerically costly. For these 
first investigations, we therefore confine ourselves to the study of correlations of the 
plaquette. 



Furthermore, we use the source method 121 to calculate the correlation function of 



the plaquettes. For this purpose, we estimate the average plaquette at the origin 
from 

Gxp I s + 'Yli{x}^iv 'nxP^uis^^\ 

W[ri\ = p ^ (15) 



JVU, exp{-^ + E 



where S is the action ([Tl|). Let P(x) denote Y^fiu P^^vix). It is straightforward to 
verify that the functional derivative of W[ri\ with respect to the source r]{x) yields 
the desired correlation function, i.e. 

c{t) = ^i.=o = {p{x)pm - {p{x)){pm . (16) 

In practice, we are interested in the correlation of the plaquette in time direction 
implying that one chooses ri{x) = ri{t). In fact, one simulates two statistical ensem- 
bles. One ensemble is generated with the inverse temperature, the other is obtained 
by setting the inverse temperature of the time slice t = to (3 + rj leaving the remain- 



ing /?- values unchanged |T^. In both ensembles, the average plaquette is obtained 
as a function of time. The correlation function (|16D is obtained by approximating 
the functional derivative in (0) by the difference of the expectation values of the 
plaquette of the two ensembles. 

We use 1140 heat bath steps to extract the average plaquette in both ensembles. A 
typical result for the correlation function as a function of time is shown in figure 
^, where the noise suppression factor is set to j = 0.5. The inverse temperature is 
f3 = 2.1 guaranteeing that the systems are in the scaling region (see left picture of 
figure |l|). One clearly observes an exponential decay of the correlation, where the 
slope of \nC(t) provides the SU(2) mass gap in units of the lattice spacing. 
It is interesting to compare the value of the mass gap, obtained with and without 
noise suppression, for /3-values corresponding to the middle and the onset of the 
scaling window, respectively. In the latter case, the contribution of "irrelevant oper- 
ators" to the mass gap should be more pronounced compared with the former case. 
The numerical results are summarized in the following table. 
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Figure 3: The correlation function (IG) as a function of time t in units of lattice 
spacings for P = 1.9 and j = 0.5. 
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The error bars indicate the uncertainty due to statistical fluctuations. They are ex- 
tracted from the fit of the correlation function C{t) to the function const. exp{— mt}. 
The string tension k sets the scale. We use k = 440 MeV. L is the extension of our 
lattice in each direction. One observes that the discrepancy between the values 
m/i/zt with and without noise suppression decreases, if the ensemble turns towards 
the continuum limit, i.e. na'^ decreases. This shows that the infiuence of the "ir- 
relevant" operators, which distinguishes our action ([l6|) from the standard Wilson 
action diminishes. Note, however, that at physical lattice sizes where the infiuence 
of the "irrelevant" operators is small, finite size effects might play a role. This would 
imply that a large number of lattice points (here we use 10^ lattice points) is neces- 
sary for high precision measurements of the SU(2) mass gap. The numerical result 
for the mass gap in the scaling region is in agreement with the results of [§. 
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4 Conclusions 



We have shown that constraining the plaquette to its average value can be un- 
derstood as composite field renormalization. The parameter j, which enters the 
action, acts as a renormalization constant absorbing the divergences arising from 
the composite nature of fields. These divergences are responsible for the small 
signal-to-noise-ratio in correlation functions of plaquettes. The action (0) with 
noise suppression differs from the standard Wilson action up to a shift in j3 only by 
"irrelevant operators", i.e. both actions coincide in the naive continuum limit. 
We have numerically studied the new action ([TT|) which implements the noise sup- 
pression on a coarse grained lattice consisting of 10"^ lattice points. We have used 
a P independent noise suppression factor j. Other choices are possible and per- 
haps more convenient depending on the type of correlation function which is under 
consideration. 

We have extracted the Creutz ratios from the numerical data with and without 
noise suppression and have established a scaling window in both cases. We have 
numerically confirmed that the action (|11|) reproduces the scaling of the standard 
Wilson action with (3 shifted to a lower value in agreement with the analytical result. 
The SU(2) mass gap has been estimated from the plaquette-plaquette correlation 
function. We have found that the contributions from the "irrelevant" operators 
to the screening mass decrease with increasing values of jS. The goal of the noise 
suppression has mainly been the reduction of the statistical fiuctuations of the "back- 
ground" , on top of which the signal exists. 

For high precision measurements of glue-ball masses, one should use correlation 
functions of operators which have a larger overlap with the glue-ball wave function 
than the plaquettes. In addition, "perfect" actions will help to extrapolate to the 
continuum limit. A generalization of the noise suppression introduced in the present 
paper to either case seems feasible. In the case of the "perfect" actions, one has 
to ensure that the noise suppression term does not spoil the correct ultra-violet 
behavior exploited by the "perfect" action technique. 
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